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I.  Introduction  and  Program  Objective 

Current  methods  of  lifetime  assessment  leave  much  to  be  desired.  Typically,  the  expected  life  of 
a  full-scale  component  exposed  to  a  complex  environment  is  based  upon  empirical  interpretations 
of  measurements  performed  on  microscopic  samples  in  controlled  laboratory  conditions. 
Extrapolation  to  the  service  component  is  accomplished  by  scaling  laws  which,  if  used  at  all,  are 
empirical;  little  or  no  attention  is  paid  to  synergistic  interactions  between  the  different  components 
of  the  real  environment.  With  the  increasingly  hostile  conditions  which  must  be  faced  in  modem 
aerospace  applications,  improvement  in  lifetime  estimation  is  mandated  by  both  cost  and  safety 
considerations. 

This  program  aims  at  improving  current  methods  of  lifetime  assessment  by  building  in  the 
characteristics  of  the  micro- mechanisms  known  to  be  responsible  for  damage  and  failure.  The 
broad  approach  entails  the  integration  and,  where  necessary,  augmentation  of  the  micro-scale 
research  results  currently  available  in  the  literature  into  a  macro-scale  model  with  predictive 
capability. 

In  more  detail,  the  program  will  develop  a  set  of  hierarchically  structured  models  at  different 
length  scales,  from  atomic  to  macroscopic,  at  each  level  taking  as  parametric  input  the  results  of 
the  model  at  the  next  smaller  scale.  In  this  way  the  known  microscopic  properties  can  be 
transported  by  systematic  procedures  to  the  unknown  macro-scale  region.  It  may  not  be  possible 
to  eliminate  empiricism  completely,  because  some  of  the  quantities  involved  cannot  yet  be 
estimated  to  the  required  degree  of  precision.  In  this  case  the  aim  will  be  at  least  to  eliminate 
functional  empiricism.  Restriction  of  empiricism  to  the  choice  of  parameters  to  be  input  to  known 
functional  forms  permits  some  confidence  in  extrapolation  procedures  and  has  the  advantage  that 
the  models  can  readily  be  updated  as  better  estimates  of  the  parameters  become  available. 

II.  Program  Organization 

The  program  has  been  organized  into  specific  tasks  and  subtasks  as  follows. 

Task  100.  Lifetimes  of  metallic  dispersed-phase  composites 

Most  service  materials  fall  into  the  category  of  dispersion-hardened  metallic  composites.  This  task 
will  consider  the  problem  of  dispersion  hardened  materials  in  general,  but  with  two  specific 
materials,  NiAl  and  MoSij/SiC  in  mind. 

Task  110.  Identification  and  modelling  of  micromechanisms 

The  purpose  of  this  task  is  to  determine  what  micromechanisms  are  operative  in  the  high- 
temperature  deformation  of  dispersion-hardened  materials.  In  the  general  case  this  will  be  done 
by  a  literature  search.  For  specific  materials,  the  micromechanisms  will  be  determined  from  the 
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experimental  program  at  NRL.  Once  identified,  each  of  these  micromechanisms  will  be  modelled, 
in  order  to  determine  what  are  the  critical  parameters  which  determine  its  effect  on  plastic  flow 
and  values  for  these  parameters.  Also  to  be  determined  is  whether  the  modelled  critical  values 
are  dependent  on  quantities  which  must  be  obtained  from  a  smaller  scale  model. 

Task  111.  Equiaxed  dispersoids 

This  task  will  consider  dispersions  of  the  type  encountered  in  NiAl-like  materials.  That  is,  the 
dispersoids  are  considered  to  be  small  compared  to  the  grain  size.  The  term  'equiaxed'  is  used 
because  the  particles  are  roughly  of  the  same  size  in  all  three  dimensions.  However,  this  is  not  a 
requirement  for  this  task.  Rather,  it  is  necessary  that  the  particles  not  be  too  large  in  the  dimension 
normal  to  the  slip  plane,  so  that  they  can  be  surmounted  with  relative  ease  by  cross-slip  and/or 
climb  without  the  generation  of  appreciable  back-stress. 

Task  112.  Anisotropic  dispersoids 

This  task  covers  the  case  of  dispersoids  which  are  elongated  in  the  direction  normal  to  the  slip 
plane.  An  example  is  SiC  fibers  in  MoSij.  In  this  case,  plastic  flow  around  the  dispersoids  takes 
place  by  a  combination  of  glide  and  climb,  but  is  a  protracted  process  during  which  large  stresses 
acting  in  opposition  to  the  applied  load  are  developed. 

Task  113.  Grain  boundary  effects 

This  task  will  examine  the  role  of  grain  boundary  processes  in  high-temperature  deformation. 

Task  120.  Macroscopic  stochastic  model  for  creep 

In  real  materials  it  is  likely  that  more  than  one  mechanism  will  be  operative,  either  in  parallel  or 
in  series.  The  information  gained  in  task  110  is  not  sufficient  to  describe  this  situation.  Once  the 
critical  parameters  for  individual  mechanisms  have  been  determined,  it  is  necessary  to  combine 
them  in  a  macroscale  stochastic  model.  This  will  be  done  by  determining  critical  stresses  and 
activation  enthalpies  as  a  function  of  local  geometry  and  using  these  values  in  a  finite-temperature 
simulation  of  creep  through  a  random  array  of  dispersoids.  Careful  attention  must  be  paid  to 
possible  interactions  between  mechanisms. 

Task  130.  Extension  to  cyclic  deformation 

The  final  step  in  task  100  is  to  extend  the  results  to  the  case  of  cyclic  deformation.  Irreversibility 
is  an  intrinsic  feature  of  the  model  in  task  120.  However,  it  is  likely  that  other,  as  yet 
unrecognized,  characteristics  of  cycled  deformation  will  have  to  be  considered. 
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Task  200.  Lifetimes  of  piezoelectric  ferroelectrics 

Failure  in  cyclic  loading  of  sensors  and  actuators  formed  from  lead  zirconate  titanate  (PZT)  is  a 
continuing  problem.  PZT  is  a  ceramic  and  therefore  differs  from  the  materials  considered  in  task 
100  in  that  plastic  deformation  is  not  involved.  This  task  will  examine,  modelling  as  necessary, 
the  operation  of  PZT  devices,  in  order  to  determine  the  factors  governing  lifetime  limitation. 

Task  300.  Reporting 

Running  concurrently  with  tasks  100  and  200,  this  task  will  inform  the  Navy  Program  Manager 
and  Contracting  Officer  of  the  technical  and  fiscal  status  of  the  program  through  R«&D  status 
reports. 
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III.  Technical  Progress 

Task  100.  Lifetimes  of  metallic  dispersed-phase  composites 

Task  110.  Identification  and  modelling  of  micromechanisms 

Three  publications  of  varying  degrees  of  relevance  to  the  present  program  are  appended  to  this 
report. 

"The  Flow  Stress  of  Potassium"  by  Duesbery  and  Basinski,  demonstrates  that  quantitative 
agreement  between  theory  and  experiment  can  be  obtained  at  the  microscopic  level  from  first 
principles  calculations  for  sophisticated  maechanical  properties  such  as  the  flow  stress. 

"Self-Repair  of  Monolayers  with  Vacancy  Damage"  by  Joos  and  Duesbery,  introduces  a  novel  and 
potentially  useful  property  of  certain  thin  coating  materials. 

"On  Hardening  by  Dispersed  Spherical  Particles"  by  Louat  and  Duesbery,  covers  the  analytic 
work  on  dispersion  hardening  which  was  detailed  in  the  last  report. 


Task  200.  Lifetimes  of  piezoelectric  ferroelectrics 

The  utility  of  PZT  as  a  'smart  material'  is  severely  limited  by  its  premature  mechanical  failure 
when  subjected  to  cyclic  loading.  The  cause  of  this  failure  forms  the  critical  question.  The  central 
problem  is  that  an  irreversible  mechanism  is  necessary  to  explain  the  experimental  observations, 
while  the  overt  physical  properties  -  that  is,  piezoelectricity,  ferroelectricity  and  ferroelasticity  - 
are  all  symmetric  to  reversal.  A  theory  will  be  presented  which  introduces  an  intrinsic 
irreversibility  and,  at  the  same  time,  explains  all  salient  experimental  observations.  Critical 
experiments  and  modelling  directions  will  be  detailed. 

The  key  experimental  observations  are  as  follows. 

★In  the  poled  state,  ferroelectric  domains  are  elongated  strongly  in  the  direction  of  the  poling  field 
and  are  separated  by  walls  with  little  curvature. 

★Under  cyclic  mechanical  loading  at  70%  of  the  monotonic  fracture  stress,  the  lifetime  is  in  excess 
of  10’  cycles.  At  stress  amplitudes  from  75  %  to  95  %  of  the  fracture  stress,  the  lifetime  decreases 
sharply  from  10*  to  10’  cycles)*.  The  effect  of  temperature  is  unknown. 

★Under  cyclic  (resonant)  electrical  loading,  the  lifetime  depends  on  temperature.  For  deformation 
above  80°  C,  failure  is  much  more  rapid  than  below  tfiis  critical  temperature.  Predeformadon 
below  80°  C  prior  to  high-temperature  deformadon  leads  to  even  more  rapid  failure.  Damage  is 
not  detectable  for  deformadon  below  80°  C,  but  is  observed  as  intergranular  cracking  at  higher 
temperatures*. 
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a.  Theory 

a.l  Domain  Wall  Motion 

There  are  two  regimes  of  domain  wall 
motion,  depending  on  the  domain 
structure  and  the  misfit  strain.  For  thin 
domains  and/or  small  misfit,  the 
incommensuracy  is  accommodated  by 
coherency  strains.  For  thick  domains 
and/or  large  misfit  strain,  the  do'^ain 
walls  are  expected  to  bind  strongly  to 
dislocations.  This  report  will  be 
concerned  only  with  the  second  regime, 
in  which  the  domain  wall  motion  is 
limited  by  dislocation  motion. 

The  dislocations  of  importance  will 
have  Burgers  vectors  and  line  directions 
lying  in  the  plane  of  the  domain  walls. 

Motion  of  the  domain  walls  normal  to 
themselves  therefore  requires  that  these  dislocations  move  by  climb.  In  general,  dislocation  climb 
is  a  process  requiring  the  emission  or  absorption  of  simple  vacancies.  For  PZT,  however,  the 
situation  is  more  complicated,  as  can  be  seen  from  Figure  1,  which  shows  the  atomic  structure 
(schematic  only)  of  a  domain  wall  dislocation  in  PZT.  The  position  of  the  wall,  which  is  the  glide 
plane  of  the  dislocation,  is  shown  in  Figure  1  as  a  horizontal  dashed  line.  It  is  clear  that  an 
upward  (downward)  i.iotion  of  the  domain  wall  requires  that  the  extra  half-plane  of  the  dislocation 
emit  (absorb)  PbOj  or  Zr(Ti)02  complexes,  which  will  require  higher  activation  energy  than  the 
simply  vacancy  process.  The  possibility  also  exists  that  the  defects  involved  in  climb  may  be 
electrically  charged,  therefore  interacting  directly  with  the  internal  electric  field. 

The  creation  of  a  climb  defect  is  a  thermally  activated  process  requiring  an  energy  U.  The 
necessary  enthalpy  for  a  unit  climb  step  of  a  domain  wall  dislocation  is 

H  =  U~xb^  (1) 


where  xb^  is  the  force  applied  to  a  length  b  of  dislocation  which  moves  through  the  same  distance 
b  as  the  climb  defect  is  created;  it  is  assumed  in  (1)  that  electrical  neutrality  is  maintained 
throughout  the  climb  process.  The  rate  of  creation  of  climb  defects  R  in  these  circumstances  is 


R 


vcxp[- 


U-xh^ 

kT 


(2) 
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where  v  is  the  frequency  of  atomic  vibration,  k  is  Boltzmann's  constant  and  T  is  the  absolute 
temperature. 


Assume  that  the  dislocations  are  spaced  a  distance  X  apart,  so  that  the  stress  is  given  by 


X  =  F— 
b 


(3) 


where  F  is  the  uniform  pressure  tending  to  move  the  wall,  given  in  the  last  report  as 

-  =  l.SOeV-)  -  LSeeV?  e— 
ft  3  fi  _ 


(4) 


in  which  f  is  the  volume  fraction  of  material  with  c-axis  aligned  along  the  poling  direction,  e  is 
the  misfit  and  o  is  the  external  stress  applied  normal  to  the  domain  wall  (positive  if  tensile). 
Substituting  (4)  and  (3)  into  (2),  and  noting  that  €=X/b,  the  velocity  of  the  wall  can  be  expressed 
as 


V  =  iviexpf - — — ] 

kT 


(5) 


where  the  positive  (negative)  signs  apply  to  tensile  (compressive)  applied  stress.  For  cyclic 
stressing,  the  net  migration  rate  (in  cm.  per  cycle)  of  the  domain  wall  is 


<v>  =  2vftexp[--^— ]sinh[  ] 

kT  kT 


(6) 


in  such  a  sense  as  to  make  the  crystal  extend  along  the  stress  axis.  If  U  is  of  order  leV,  equation 
(6)  predicts  a  wall  drift  of  40yxm  in  10’  cycles.  It  is  expected  that  U  will  be  substantially  larger, 
2eV  or  more,  and  in  this  case  the  wall  motion  will  be  negligible  (10“  cycles  needed  for  boundary 
drift  of  Inm). 
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b.  Plans 

*  Continued  modelling  of  the  mechanisms  of  ferroelastic  domain  wall  motion  will  be  performed, 
with  the  aim  of  determining  critical  parameters  which  can  be  derived  from  or  compared  with 
experiment.  Specific  questions  to  be  answered  are: 

1.  Is  the  wall  motion  governed  by  the  two-dimensional  analog  of  kink  pair  nucleation  and 
propagation,  or  by  diffusive  dragging  of  wall  defects  and  dislocations? 

2.  How  are  these  mechanisms  influenced  by  an  electric  field  of  sufficient  magnitude  to  stabilize 
the  poled  state? 

3.  How  do  the  lifetimes  of  poled  and  unpoled  specimens  compare? 
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Abstract — Experimental  results  on  the  flow  stress  of  potassium  in  the  temperature  range  of  1-30  K  are 
compared  with  atomistic  computer  simulation  of  rigid  dislocation  motion  and  of  dislocation  kink  pair 
generation  and  migration.  Quantitative  agreement  is  demonstrated.  Implications  for  the  rate  theory  and 
the  continuum  elastic  theory  of  kink  pairs  are  discussed. 


1.  INTRODUCTION 

It  is  now  accepted  that  the  large,  strongly  tempera¬ 
ture  and  orientation  dependent  flow  stress  exhibited 
by  even  the  simplest  body-centred  cubic  (b.c.c.)  crys¬ 
talline  materials  has  its  origin  in  a  symmetry-based 
non-planar  structure  of  the  rate-controlling  screw 
dislocation.  Suggested  first  by  Hirsch  [I],  this  was 
established  by  extensive  atomistic  simulations  which 
have  been  reviewed  in  depth,  most  recently  by  Vitek 
[2]  and  Duesbery  (3]. 

In  the  subsequent  shift  of  focus  of  atomistic  simu¬ 
lation  to  more  complex  materials,  it  has  been  forgot¬ 
ten  that  some  basic  but  readily  addressable  questions 
remain  to  be  answered.  To  illustrate  this  more  clearly, 
brief  recourse  to  the  history  of  the  subject  is  in  order. 
During  the  early  1970s,  before  the  advent  of  the 
embedded  atom  method  and  its  adjuncts  [4, 5],  the 
almost  universal  use  of  pair  potentials  led  to  two 
distinct  approaches.  The  first  of  these,  acknowledging 
that  the  refractory  b.c.c.  metals  cannot  be  represented 
with  any  confidence  by  pair  potentials,  concentrated 
instead  on  investigating  the  range  of  dislocation  core 
properties  which  could  be  supported  by  a  b.c.c.  lattice 
bound  by  empirical  pair  potentials.  This  avenue  of 
research  was  used  with  great  success  by  Vitek  and 
co-workers  [2, 3].  They  demonstrated,  among  other 
things,  the  universality  of  the  large  Peierls  binding 
energy  of  the  screw  dislocation.  The  second  ap¬ 
proach,  followed  by  Basinski  and  colleagues  [2,3], 
recognised  that  pair  potentials  can  be  used  with  a 
high  degree  of  justification  for  alkali  metals  and 
therefore  that  a  combined  theoretical  and  exper¬ 
imental  attack  on  the  properties  of  potassium, 
the  only  alkali  metal  which  does  not  experience  a 
low-temperature  martensitic  transformation,  might 
permit  quantitative  agreement  to  be  realised.  It  is 
with  this  second  approach  that  this  paper  will  be 
concerned. 

The  work  available  on  potassium  is  quite  extensive. 
The  core  structure  and  mobility  of  rigid  screw  dis¬ 


locations  in  the  presence  of  shear  stresses  were  inves¬ 
tigated  theoretically  by  Basinski  et  at.  [6).  This  work 
was  extended  to  encompass  general  applied  stresses 
by  Basinski  and  Duesbery  [7]  and  by  Duesbery  [8]. 
The  question  of  non-rigid  screw  dislocations  was 
addressed  by  Duesbery  [9],  who  investigated  the 
structure  and  mobility  of  kink  pairs  in  potassium. 
Meanwhile,  the  temperature  and  orientation  depen¬ 
dence  of  the  flow  stress  in  potassium  for  temperatures 
between  I  and  30  K  were  measured  by  Basinski  et  at. 
[10].  The  results  of  this  work  were  never  fully 
integrated  with  the  original  purpose  of  quantifying 
computer  modelling  of  the  dislocation  core.  The 
purpose  of  this  paper  to  correct  this  omission.  It  will 
be  shown  that  when  reliable  interatomic  forces  can  be 
used,  a  theoretical  treatment  of  the  screw  dislocation 
core  is  sufficient  to  give  a  quantitative  description  of 
the  orientation  and  temperature  dependence  of  plas¬ 
ticity  in  a  b.c.c.  metal.  Some  associated,  fundamental 
questions  will  also  be  discussed. 

2.  THE  RIGID  SCREW  DISLOCATION 

The  initial  work  on  the  mobility  of  screw  dis¬ 
locations  in  potassium  [6]  was  performed  for  rigid, 
straight  dislocations  at  zero  temperature,  a  restriction 
common  at  the  time  to  all  core  simulations  and 
imposed  by  computational  limitations.  Using  pure 
shear  strains  to  apply  a  force  to  the  core,  the  Peierls 
stress  was  found  to  lie  between  1.05  x  10“’ /i  (/t  is 
the  <lll>{fik/}  shear  modulus)  for  a  strain  applied 
parallel  to  the  easy  direction  on  {211}  planes  and 
2.2  X  10“’ /i,  for  strains  applied  on  {110}  planes.  It 
was  found  later  [6, 7]  that  the  core  in  b.c.c.  metals  is 
particularly  sensitive  to  the  details  of  the  applied 
stress  tensor;  calculations  for  forces  applied  with 
uniaxial  tensile  loading  revised  the  limits  on  the 
calculated  Peierls  stress  to  between  6.3  x  10“’ /i  (for 
<110>  loading)  and  8.5  x  10“  V  (for  loading  parallel 
to  an  axis  in  the  centre  of  the  standard  stereographic 
triangle). 
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The  experimental  work  [lOJ,  however,  showed 
Peierls  stresses  varying  only  from  2.0  x  10  (for 
<II0>  loads)  to  33xl0'';i  (for  centre-triangle 
loads),  even  when  extrapolated  to  0  K.  This  dis- 
erepaney.  a  factor  of  3.  is  well  outside  either  the 
experimental  error  or  the  theoretical  uncertainty  and 
calls  into  question  the  validity  of  the  mechanism 
assumed  to  be  operative  in  the  computer  simulations, 
i.c.  the  surmounting  of  the  Peierls  barrier  by  a  riipd 
dislocation.  At  the  same  time,  the  experimental 
measurements  of  the  orientation  dependence  of  the 
flow  stress  provided  corroborative  evidence  support¬ 
ing  the  simulations.  It  was  found  that  the  flow  stress 
r.  a  function  of  the  temperature  T and  the  orientation 
G,  could  be  represented  in  the  separable  form 

t(r.e)=/(r)  g(0).  (D 

This  means  that  the  orientation  dependence  of  the 
flow  stre.ss  is  an  explicit  function  only  of  the  magni¬ 
tude  of  the  stress  and  not  of  the  temperature.  In  other 
words,  the  ratio  of  the  flow  stresses  for  two  different 
orientations  at  the  same  temperature  is  independent 
of  temperature.  With  the  temperature  factored  out  in 
this  way,  it  was  found  [10]  that  the  theoretical  and 
experimental  orientation  depcnrences  of  the  flow 
stress  in  potassium  were  in  close  agreement  (to  within 
better  than  ±5%).  It  was  shown  later  [3]  that  this 
same  reduction  (I)  could  be  applied  to  the  exper¬ 
imental  measurements  of  the  orientation  dependence 
of  the  flow  stress  in  refractory  b.c.c.  metals  with 
considerable  success  (standard  deviations  reduced  in 
all  cases  but  one  to  better  than  10%  from  up  to  30%). 

The  picture  presented  at  this  stage  was  of  a  mod¬ 
elling  procedure  which  agreed  well  with  experiment  in 
one  respect,  the  orientation  dependence  of  the  Peierls 
stress,  but  which  was  sadly  lacking  in  its  prediction 
of  the  magnitude  and  temperature  dependence  of  the 
Peierls  stress. 

3.  KINK  PAIRS  AND  THERMAL  ACTIVATION 

The  generation  and  propagation  of  kink  pairs  in 
potassium  were  modelled  by  Ducsbery  [9],  using  a 
model  containing  10'  atoms  arranged  in  a  block  of 
dimensions  306  x  306  (6  is  the  magnitude  of  the 
Burgers  vector)  normal  to  the  dislocation  line  and 
IS06  parallel  to  the  dislocation  line.  In  what  was 
strictly  a  zero  temperature  calculation,  kink  pairs  of 
varying  separation  were  equilibriated  with  respect  to 
an  applied  stress  of  such  a  sense  as  to  expand  the 
pairs.  From  this  stress-separation  data  it  is  possible 
to  extract  the  thermodynamic  parameters  which 
appear  in  the  standard  Arrhenius  expression  for  the 
strain  rate  i 

-  (Ht»  -  T*<1 

<■  =foC''''*'^  =  tue  "  (2) 

in  which  k  is  Boltzmann's  constant.  T  is  the  tempera¬ 
ture  and  if  a  pre-exponcntial  factor  to  be  discussed 
below.  As  indicated  in  (2),  the  activation  enthalpy  // 
for  a  kink  pair  in  the  presence  of  an  applied  stress  r 
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Fig.  I .  The  equilibrium  stress  and  configurational  energy  of 
a  kink  pair  complex  in  potassium  as  a  function  of  kink 
separation. 


consists  of  two  parts.  The  first  is  the  configurational 
energy  £(t)  of  the  complex.  This  term  is  reduced  by 
the  external  work  done,  which  is  the  product  of  the 
applied  stress  r,  the  area  A  swept  by  the  activation 
process  and  the  displacement  6  through  which  the 
area  is  translated. 

The  stress-separation  data  obtained  earlier  [9]  are 
reproduced  in  Fig.  1.  The  scatter  in  the  data  points 
is  larger  than  is  usual  for  theoretical  work.  This  is 
indicative  of  the  diflrculty  of  the  calculations.  The 
equilibrium  states  are  unstable,  so  that  the  isolation 
of  their  location  requires  the  application  of  a  trial 
bias  stress  together  with  a  superimposed  oscillating 
stress  of  steadily  decreasing  amplitude.  The  curve 
drawn  through  the  data  points  is  a  least  squares 
spline  fit.  For  small  jcparaiions.  the  stress  is  nearly 
constant.  In  this  regime,  as  discussed  earlier  [9],  the 
kink  pair  is  simply  a  bulge  in  the  dislocation  line,  with 
height  less  than  the  distance  between  two  Peierls 
valleys.  The  height  and  width  of  the  bulge  increase  in 
constant  proportion  until  the  adjacent  Peierls  mini¬ 
mum  is  reached.  The  kinks  then  separate  according 
to  classical  expectations  with  consequent  decrease  in 
stress.  The  energy  curve  in  Fig.  I  was  obtained  by 
integration  of  the  stress-separation  curve;  the  energy 
values  are  slightly  higher  than  those  reported  pre¬ 
viously  [9],  but  arc  more  accurate,  in  good  agreement 
with  the  asymptotic  limit  of  twice  the  energy  of  an 
isolated  kink,  calculated  from  a  model  containing 
only  a  single  kink  to  be  0.038  eV  [9). 

The  information  in  Fig.  I  can  be  reprocessed  to 
give  the  activation  enthalpy  as  a  function  of  stress. 
The  energy  £(t)  is  taken  directly  from  the  relevant 
curve  and  the  swept  area  A  is  estimated  by  multi¬ 
plying  the  kink  pair  separation  by  the  kink  height  a. 
This  latter  procedure  is  not  accurate  for  large  stresses, 
when  the  kink  pair  height  does  not  reach  its  full 
extent,  but  this  does  not  affect  the  arguments  whidi 
follow.  The  result  is  shown  in  Fig.  2.  The  enthalpy 
drops  from  a  high  of  0.076  eV  at  zero  stress  [when  the 
work  term  in  (2)  is  identically  zero]  to  0  for  a  kink 
pair  separation  of  about  126.  corresponding  to  a 
stress  of  4  X  I0~'  p  (at  this  stress  level  the  use  of  the 
full  kink  height  in  the  area  calculation  is  accurate). 
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Fig.  2.  Logarithmic  plot  of  stress  against  kink  pair  separ¬ 
ation. 

This  stress  in  only  half  that  required  to  move  a  rigid 
dislocation  of  0  K  and  leads  to  the  conclusion, 
contrary  to  conventional  expectation,  that  the  motion 
of  rigid  dislocations  may  be  superseded  even  at  the 
lowest  temperatures  by  kink  pair  activation. 

The  experimental  work  [10]  indicates  that  the 
activation  enthalpy  is  closely  proportional  to  the 
temperature  for  temperatures  below  20  K.,  with  a 
slope  of  IQk.  Using  this  single  piece  of  empirical 
information  to  fix  a  value  for  the  ratio  e/e'o  in  (2),  the 
theoretical  temperature  dependence  of  the  flow  stress 
T  can  be  estimated  from  the  data  in  Fig.  2.  The  results 
are  shown  in  Fig.  3.  together  with  the  experimental 
measurements  from  [10],  each  normalised  to  the  flow 
stress  at  4.2  K.  Agreement  is  satisfactory  up  to  20  K. 
Above  this  temperature  the  comparison  is  invalid, 
because  the  experimental  proportionality  between 
H  and  T  breaks  down.  There  remains  a  small  dis¬ 
crepancy  between  the  absolute  values  of  the  exper¬ 
imental  and  theoretical  flow  stresses.  For  an 
orientation  corresponding  to  the  constraint  used  in 
the  kink  modelling  work  [9]  (a  pure  shear  strain 
acting  on  <lll){110}),  the  experimental  values  at 
4.2  K  range  from  2.8  to  3.3  x  10  "’/r.  The  theoretical 
estimate  is  3.9  x  lO"’/!-  about  20%  too  high.  Never¬ 
theless,  the  overall  agreement  demonstrates  that  com¬ 
puter  simulation  methods  can  give  quantitative 
results  for  structures  in  which  the  interatomic  forces 
can  be  reliably  approximated. 

An  obvious  question  arises  as  to  why  the  rigid 
dislocation  core  calculations  were  found  adequate 


Fig.  3.  The  normalised  theoretical  and  experimental  tem¬ 
perature  dependence  of  the  flow  stress  in  potassium. 


for  calculation  of  the  orientation  dependence  of  the 
flow  stress.  In  the  absence  of  definitive  kink  pair 
calculations  as  a  function  of  general  applied  stresses, 
this  remains  a  matter  of  informed  conjecture.  The 
mobility  of  the  rigid  dislocation  depends  strongly 
on  the  shape  of  the  core,  which  depends  in  turn  on 
the  detailed  nature  of  the  applied  stress.  In  other 
words,  the  magnitude  of  the  Peierls  potential  depends 
critically  on  the  external  stress.  Since  the  energy  of 
the  kink  pair  complex  also  depends  on  the  Peierls 
potential,  it  is  not  unreasonable  to  suppose  that  the 
rigid  dislocation  orientation  dependence  carries 
through  to  kink  pair  activation. 

4.  CRITICAL  QUA.VnriES 

The  calculations  reported  in  the  previous  section 
are  self-consistent,  independent  of  both  elasticity 
theory  and  experiment,  except  for  the  use  of  the 
experimentally  observed  value  and  constancy  of  the 
exponent  HjkT.  The  results  can  therefore  be  used 
with  some  degree  of  confidence  to  examine  some 
associated,  critical  questions  and  quantities. 

First  to  be  considered  will  be  the  pre-exponential 
factor  4,  which  can  be  decomposed  according  to 
standard  rate  theory  into 

4  =  jAoAv(2)  (3) 

where  p  is  the  dislocation  density,  a  is  the  kink  pair 
height,  v(,i)  is  (he  attack  frequency  at  (he  critical 
activation  length  ).  and  A  is  the  kink  mean  free  path. 
As  noted  above,  experiment  indicates  a  value  for 
ln(£7£o)  of  20.  The  question  is  whether  or  not  this 
number  is  consistent  with  reasonable  magnitudes  for 
the  physical  quantities  involved  The  experiments  [10] 
were  performed  for  a  tensile  strain  rate  of 
3.4  X  10'*  s  '.  Taking  half  this  value  as  the  shear 
strain  rate,  a  value  of  8.25  x  10*  is  obtained  for  4- 
Using  a  value  of  0.3  THz  for  the  frequency  of  trans¬ 
verse  phonons  with  wave  vector  parallel  to  <  1 1 1  >  and 
wavelength  of  15b  (a  %  12.5h)  in  potassium  [1 1]  and 
setting  as:b  =  0.475  nm,  consistency  requires  the 
product  pA  to  have  a  value  of  1.44  x  10*  m"'.  With 
a  reasonable  mobile  dislocation  density  at  yield 
of  10"  m*^  the  kink  mean  free  path  must  re  300h, 
an  acceptable  value.  Participation  of  vibrational 
modes  with  larger  wavelength  (the  0.3  THz  mode, 
with  Av  a;  2.6  kT at  4  K,  is  partially  quenched)  would 
reduce  this  predicted  mean  free  path. 

The  second  critical  quantity  to  be  considered  con¬ 
cerns  the  clastic  interaction  of  kinks,  the  continuum 
elastic  theory  of  which  indicates  that  the  two  mem¬ 
bers  of  a  kink  pair  of  height  a  and  separation  /  attract 
each  other  with  a  force  given  by  [12] 

where  1^(6)  is  the  anisotropic  elastic  energy  factor  of 
a  dislocation  line  of  character  angle  0  (i.e.  such  that 
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the  energy  per  unit  length  of  the  dislocation  line  is 
W(d)b^  ln(^/p),  where  R  and  p  are  outer  and  inner 
radii).  The  expression  (4)  is  valid  provided  that  the 
kink  separation  is  much  larger  than  the  kink  width. 
While  this  is  not  true  in  the  present  case,  since  the 
maximum  pair  separation  is  only  about  four  times  the 
kink  width,  it  might  be  expected  that  (4)  would  be  a 
fair  approximation.  A  more  important  drawback  of 
(4)  is  the  assumption  that  standard  linear  elasticity 
can  be  applied  to  dislocation  segments  with  length  of 
the  same  order  as  the  inner  cutoff  radius.  For  b.c.c. 
metals,  in  which  the  screw  dislocation  core  is  known 
to  be  highly  structured  over  distances  of  this  order 
and  in  which  the  kinks  pass  over  a  Peierls  energy 
barrier,  this  assumption  is  in  particular  doubt. 

There  is  sufficient  information  in  the  data  outlined 
above  to  test  the  applicability  of  (4).  In  the  saddle 
point  state,  the  force  (4)  can  be  equated  to  the  force 
lab  exerted  on  a  kink  by  the  applied  stress.  The 
resulting  formula  relates  the  stress  directly  to  the  kink 
pair  separation  and  can  be  tested  using  the  data  from 
Fig.  I .  The  results  are  shown  in  Fig.  4.  The  straight 
line,  drawn  to  have  a  slope  of  —2  and  positioned  by 
eye,  shows  that  the  inverse  square  law  of  (4)  is  obeyed 
well  for  stresses  less  than  4.5  x  10"’  p.  It  was  re¬ 
ported  incorrectly  [9],  as  pointed  out  by  Kirchner 
[13],  that  the  intercept  of  this  line  was  also  compatible 
with  (4).  In  fact,  this  is  not  the  case  at  all.  The 
intercept  of  the  line  at  In  t  =  0  occurs  at  /  =  22.9i). 
Since  the  kinks  are  very  wide,  it  is  reasonable  to 
assign  to  them  the  stiffness  of  a  screw  dislocation, 
that  is,  [iV(0)  +  IV'iO)]  =  5.3  GPa.  Using  these  quan¬ 
tities  with  the  value  of  /r  —  1 .3  GPa  in  (4),  an  effective 
kink  height  of  a  =  2.6b  is  obtained. 

This  source  of  thi;  discrepancy  coulo  lie  in  an 
insufficiency  of  the  expression  (4)  or  in  the  assump¬ 
tion  that  the  force  on  a  kink  is  xba.  The  latter  is 
unlikely,  since  it  is  based  on  thermodynamic  argu¬ 
ments.  It  must  be  concluded  that  the  elastic  force 
acting  between  two  kinks  in  the  b.c.c.  structure  is  of 
the  same  functional  form  as  (4),  but  of  substantially 
larger  magnitude.  Interest.ng  in  this  respect  is  a  recent 
suggestion  by  Seeger  [14].  In  analyses  of  the  tempera¬ 
ture  dependence  of  the  flow  stress  in  refractory  metals 
in  terms  of  the  kink  pair  mechanism,  using  (4),  values 


Fig.  4.  The  activation  enthalpy  for  kink  pairs  in  potassium 
as  a  function  of  applied  stress. 


of  a  =  1.5b  to  1.8h  were  obtained  for  the  kink  height. 
From  this  deduction  it  was  concluded  that  the  pre¬ 
dominant  kink  in  b.c.c.  metals  was  the  doub/e  kink, 
extending  over  two  Peierls  hills.  Calculations  using 
both  pair  and  embedded  atom  potentials  [3]  do  not 
support  this  conclusion.  In  the  light  of  the  present 
results,  it  seems  more  likely  that  the  deduced  extent 
of  the  double  kink  is  a  manifestation  of  the  break¬ 
down  of  the  elastic  expression  (4). 

5.  SUMMARY  AND  CONCLUSIO.NS 

The  results  described  above  show  that  quantitative 
predictions  can  be  made  with  the  use  of  atomistic 
computer  modelling  methods  when  the  interatomic 
forces  are  known  with  some  precision.  It  would  be  an 
interesting  exercise  to  determine  whether  the  modem 
embedded  atom  potentials  can  provide  a  similarly 
satisfactory  description  of  the  refractory  b.c.c.  and 
other  transition  metals.  In  view  of  the  results  pre¬ 
sented  above,  this  might  be  considered  a  critical 
test  of  the  reliability  of  the  potentials  for  use  in 
calculation  of  other  properties. 

The  low-temperature  deformation  of  potassium 
has  been  shown  to  be  governed  by  the  generation  of 
fully  formed  kink  pairs.  The  mutual  interaction  of  the 
kinks  which  form  the  pair  was  found  to  have  the  same 
functional  dependence  on  separation  as  is  predicted 
by  classical  elasticity,  but  a  magnitude  larger  by  more 
than  a  factor  of  2.  The  discrepancy  is  tentatively 
attributed  to  one  or  more  of  several  factors,  including 
the  application  of  the  line  tension  model  of  a  dis¬ 
location  to  kinks  with  a  Peierls  energy  and  to  the 
atomic  dimensions  and  commensuracy  of  the  kink 
height  and  the  dislocation  core  structure.  Further 
work  is  needed  to  resolve  this  problem.  Particularly 
illuminating  would  be  the  properties  of  a  system 
controlled  by  incomplete  kink  pair  generation  (i.e.  in 
the  high-stress  regime  when  the  height  of  the  critical 
kink  complex  is  less  than  a  repeat  spacing)  and  of 
a  system  controlled  by  kink  migration  rather  than 
generation.  Both  types  of  system  are  likely  to  be 
encountered  with  refractory  b.c.c.  metals  and  with 
semiconductors  [3]. 

Most  computer  simulation  work,  to  this  day,  is 
confined  to  consideration  of  straight,  rigid  dis¬ 
locations.  This  constraint,  no  longer  imposed  by 
computational  capability,  is  difficult  to  justify  in  the 
light  of  modem  knowledge  that,  with  few  exceptions, 
dislocations  move  by  kink  pair  activation  and 
migration.  There  is  no  doubt  that  the  study  of  rigid 
dislocations  can  provide  much  useful  information  in 
return  for  limited  effort.  For  example,  the  orientation 
dependence  of  dislocation  mobility  can  be  assessed 
with  some  confidence.  However,  bridging  the  gap 
between  the  atomistic  and  macroscopic  scales  by 
means  of  linear  elasticity  can  be  a  dangerous  pro¬ 
cedure,  as  shown  above.  Larger  scale,  three-dimen¬ 
sional  atomistic  calculations  are  needed  to  establish 
a  reliable  continuum  interface. 
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Lennard'Joncs  monolayers  are  shown  by  molecular  dynamics  simulation  to  exhibit  a  self-healing  be¬ 
havior  which  involves  a  rapid  condensation  of  vacancies  into  dislocation  dipoles.  An  associated  oolUpae 
of  the  shear  modulus  similar  to  the  Kosterlitz-Thouless  dipole  dissociation  occurs  for  high  vacancy  con¬ 
centrations.  By  contrast  monolayers  having  nearest  neighbor  piecewise  linear  interactions  are  found  to 
cavitate  by  agglomeration  of  vacancies  into  clusters. 

PACS  numbers:  62.20.Dc.6l.72.Ji.68.SS.Lii 


Vacancy  damage  is  a  common  occurrence  in  mono- 
layers  and  thin  Alms,  more  so  than  in  bulk  solids  because 
thermodynamic  equilibrium  in  the  former  involves  ex¬ 
change  with  the  external  environment  as  well  as  entropic 
considerations.  Vacancies  are  particularly  abundant  in 
physisorbed  monolayers,  in  which  the  concentration  can 
be  as  high  as  10%  before  the  melting  point  is  reached  [ll. 
This  Letter  shows  that,  although  vacancies  are  local  de¬ 
fects,  a  concentration  of  just  a  few  percent  can  have  pro¬ 
found  consequences  on  the  macroscopic  physical  proper¬ 
ties  of  a  monolayer.  Specifically  we  demonstrate,  by 
molecular  dynamics  simulation,  a  unique  self-repair  me¬ 
chanism  which  leads  to  a  shear  modulus  collapse  similar 
to  the  melting  scenario  proposed  in  the  Kosterlitz- 
Thouless-Halperin-Nelson-Young  (KTHNY)  l2l  theory 
of  two-dimensional  melting.  In  this  new  mechanism,  va¬ 
cancies  condense  rapidly  into  dislocation  dipoles,  repair¬ 
ing  the  lattice;  the  dipoles  then  dissociate  inducing  a 
quasifluid  state.  These  effects  occur  for  temperatures  just 
above  the  self-diffusion  temperature. 

Monolayers  can  appear  in  many  forms  with  various  in¬ 
teractions.  We  focus  on  monolayers  formed  of  atoms  in¬ 
teracting  with  spherical  potentials.  Two  potentials  are 
considered,  one  long  range,  the  Lennard-Jones  potential 
(UP),  and  the  other  short  range,  the  nearest  neighbor 
piecewise  parabolic  potential  (PWPP).  The  effects  dis¬ 
cussed  above  are  seen  in  the  UP  monolayer.  The  other 
system,  the  PWPP  monolayer,  exhibits  very  different  be¬ 
havior.  It  does  not  heal  but  develops  voids. 

The  UP  energy  between  two  atoms  at  a  distance  r 
from  each  other  is  given  by 


Ku(r)  -4< 


12 

6' 

_ 

r 

r 

(1) 


The  PWPP  energy  between  two  atoms  at  a  distance  r 
from  each  other  is  given  by 

rfr  —  </o)  *  “  mv  *,  r^da+w, 

>'pwp(/’)'"j  —  r(r  — </o~2w)*,  <r  :Ssdti'^lw , 

0,  r></o+2w. 

(2) 


w  is  taken  to  be  O.ISc/o  as  in  Ref.  [3).  For  this  value  of  w 
the  dislocation  core  has  only  one  broken  bond.  Larger 
values  of  w  lead  to  more  broken  bonds  in  the  core  [4). 

Monolayers  formed  of  particles  interacting  with  these 
potentials  condense  into  triangular  lattices.  Both  types  of 
monolayers  have  been  extensively  studied  although  the 
motivations  of  the  studies  were  different.  The  PWPP  has 
been  used  mostly  in  simulations  of  the  plasticity  of  solids 
(3,4]  whereas  the  UP  was  used  extensively  in  attempts  to 
understand  the  nature  of  the  melting  transition  in  two  di¬ 
mensions  (for  reviews  see  Refs.  [5]  and  (61).  Recently 
the  melting  behavior  of  the  PWPP  was  also  considered 
17J. 

An  expression  for  the  elastic  constants,  which  is  valid 
at  finite  temperatures,  can  be  derived  by  noticing  that  the 
elastic  constants  are  the  second  derivatives  of  the  Helm¬ 
holtz  free  energy,  O  “  —kTlnZ,  with  respect  to  the  vari¬ 
ous  Lagrange  strains  tjij  (8]: 

The  classical  partition  function  Z  for  the  crystal  is  a  rela¬ 
tively  simple  sum  over  all  configurational  energies. 

For  an  isotropic  solid  such  as  the  two-dimensional  tri¬ 
angular  lattice,  there  are  only  two  independent  elastic 
constants.  Traditionally  the  following  two  are  chosen  to 
J  characterize  the  system: 


d*n 


(3) 
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where  the  quantities  Ax  and  Ay  represent  the  x  and  y  r 
coordinates  of  (x**— x'O,  r  the  modulus  of  that  vector, 
and  A  the  area  of  the  system.  The  terms  in  braces  are  re* 
ferred  to  as  the  "fluctuation  terms”  while  the  following 
terms  in  each  expression  constitute  the  “Born  terms.” 
The  last  term  in  Eq.  (5)  is  the  kinetic  energy  contribu¬ 
tion. 

For  canonical  ensembles  the  above  equations  ((4)  and 
(S)]  are  equivalent  to  the  ones  given  by  Ray,  Moody,  and 
Rahman  [9].  Whether  the  molecular  dynamics  simula¬ 
tion  is  done  at  constant  pressure  or  constant  volume,  the 
sums  on  the  particle  momenta  can  be  factored  out  using 
the  equipartition  theorem.  This  is  the  case  provided  the 
simulation  maintains  a  canonical  ensemble. 

To  simulate  a  physically  realistic  system  we  have  per¬ 
formed  simulations  at  constant  pressure  and  varied  the 
temperature.  A  slightly  compressed  monolayer  was 
chosen.  The  two  potentials  have  been  chosen  to  have 
the  same  depth  <6"*irw^)  and  the  same  minimum  (a 

The  zero  temperature  elastic  constants  for  a  perfect 
lattice  at  a  pressure  of  0.0053x  with  the  LIP  interaction 
are  X*0.836r  and  /i*’0.825r  whereas  with  the  PWPP 
they  are  X  •■0.440ic  and  n  “0.429x.  The  UP  system  is  a 
stifier  solid. 

Figure  1  shows  the  variation,  as  a  function  of  vacancy 
concentration  Ch,  of  the  shear  modulus  (/i)  and  the  bulk 
modulus  (B“X+p),  normalized  to  their  respective  per¬ 
fect  lattice  values  no  and  Bo.  Vacancy  concentrations 
from  0  to  3.8  at.%  were  considered  and  the  calculations 
were  done  at  a  low  temperature,  O.OOlxdi/k,  where  k  is 
‘the  Boltzmann  factor.  (N.B.  The  melting  temperature 


vaconcy  concentration  (at.  /J) 

FIG.  1.  Variation  of  the  shear  modulus  y  (squares)  and  the 
bulk  modulus  (circles),  both  normalized  to  the  zero 

vacancy  values  yo  and  Bo,  respectively,  as  a  function  of  vacancy 
concentration  for  the  UP  (full  symbols)  and  PWPP  (open  sym¬ 
bols)  systems  at  the  temperature  0.002icd3/ka. 


for  both  systems  is  about  0.01 1  xdi/k.) 

The  range  of  the  potential  seems  to  have  little  effect  on 
the  rate  of  decrease  of  the  shear  modulus  with  increasing 
vacancy  concentration.  The  shear  modulus  loses  a  fifth  of 
its  value  for  a  4  at.%  concentration  of  vacancies.  This 
corresponds  to  "PoV/ro"  “Sc,.  The  bulk 

modulus  on  the  other  hand  is  less  sensitive  to  vacancy  de¬ 
fects;  the  rate  of  decrease  is  considerably  smaller  (AB/Bo 
»  — l.2Sc,).  As  Fig.  I  shows,  there  is  a  scatter  in  the 
data  with  vacancies  present.  For  each  concentration 
there  are  many  possible  configurations  of  vacancies.  Fig¬ 
ure  1  shows  one  configuration  for  each  studied  concentra¬ 
tion. 

As  the  temperature  is  raised,  the  vacancies  start  to 
move.  In  both  systems  the  self-diffusion  temperature  for 
vacancies  is  approximately  0.0055  ±0.000Sx<fJ/A:.  The 
two  systems,  however,  show  rather  different  behavior. 

In  the  UP  system,  the  vacancies,  when  in  sufficient 
number,  condense  into  a  dislocation  dipole.  The  sample 
in  our  simulations  is  a  rhombus  with  28  atoms  on  the 
side.  Periodic  boundary  conditions  are  used.  For  fewer 
than  15  vacancies  (or  c,  <  1.9%)  the  vacancies  diffu 
into  small  vacancy  clusters.  For  1 5  vacancies  or  more  (SP 
1.9%)  the  vacancies  are  removed  by  condensing  into 
a  “large”  vacancy  dislocation  dipole  (separation  equal  to 
or  larger  than  15  atomic  rows)  [see  Fig.  2(b)]  [10,11], 
The  periodic  boundary  conditions  introduce  a  periodicity 
in  the  behavior  of  our  system  [see  Fig.  3(a)].  For  more 
than  28  vacancies  (or  c,  >  3.6%)  there  are  more  vacan¬ 
cies  than  would  fit  into  one  dipole,  so  28  vacancies  are  an¬ 
nealed  out  and  the  remainder  form  vacancy  clusters  un¬ 
less  their  numbers  exceed  the  minimum  required  to  form 
a  second  “large”  dipole.  When  a  dipole  is  present  in  the 
sample,  the  shear  modulus  collapses  [see  Fig.  3(a)].  It 
takes  about  a  time  of  300(X)ro  to  form  a  “large”  dipole, 
i.e.,  about  20  ns  if  parameters  appropriate  to  Xe  are  used 
(/q  is  the  unit  of  time  Vm/x).  This  is  an  upper  limit  ob¬ 
tained  with  vacancies  uniformly  placed  on  the  sample. 
With  other  configurations  the  times  were  shorter.  Modu¬ 
lations  in  a  substrate  holding  field  could,  however,  make 
these  times  much  longer,  leading  to  an  increase  in  the 
concentration  of  vacancies  (e.g.,  Xe  monolayers  adsorbed 
on  graphite  [|]). 

The  above  results  are  in  accord  with  our  findings  on  the 
properties  of  dislocation  dipoles  [10-12].  In  these  refer¬ 
ences  it  is  shown  that  dipoles  containing  fewer  than  15 
vacancies  behave  more  like  composite  pinned  defects  ths” 
elastic  dipoles.  They  are  the  favored  defects  at  very  Ic  ^ 
vacancy  concentrations.  On  the  other  hand  a  logarithmic 
dependence  of  the  energy  on  size  makes  large  dipoles 
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FIG.  2.  (a)  Typical  low  temperature  configuration  of  atoms 
for  either  the  LJP  or  PWPP  system  for  a  vacancy  concentration 
of  r«~2.S  at.%.  (b)  Same  configuration  after  relaxation  with 
the  LJP  interaction  at  the  temperature  T“0.006Kdi/kB. 
(c),(d)  Same  as  (b)  but  with  the  PWPP  interaction:  (c)  is  an 
intermediate  configuration  and  (d)  the  final  configuration. 


favorable  at  higher  concentrations.  These  move  nearly 
freely  in  the  direction  of  their  Burgers  vector  over  a  wide 
range  of  angles  (45"  to  135")  [12]  and  are  easily 
screened  out  by  similar  dipoles  (in  our  case  in  the  neigh* 
boring  cells  of  our  periodic  structure).  When  present 


vacancy  eoncenfration  (at.  ^ 

FIG.  3.  Variation  of  the  shear  modulus  as  a  function  of  va* 
cancy  concentration  for  temperatures  of  0.006  and  0.009nf|/J[a 
for  (a)  the  LJP  system  (at  a  pressure  P~0.0087r)  and  (b)  the 
PWPP  system  (P-0.0053)r). 


they  lead  to  the  loss  of  shear  strength  and  translational 
order. 

When  particles  interact  with  the  short>range  PWPP, 
the  vacancies  form  clusters  which  eventually  combine  into 
large  voids  [see  Figs.  2(c)  and  2(d)].  This  typically  takes 
much  longer  than  the  condensation  in  the  UP  system  be* 
cause  of  the  short  range  of  the  interaction  (I  to  2  orders 
of  magnitude  longer).  The  PWPP  is  nearest  neighbor,  so 
as  soon  as  a  triangle  of  vacancies  is  formed,  additional 
vacancies  in  the  cluster  do  not  add  to  the  energy.  In 
ciTect,  one  has  a  free  surface,  and  the  vacancies  anneal 
out  but  remain  within  the  system.  The  temperature- 
induced  diffusion  smooths  out  the  elastic  constants  even 
before  the  voids  are  formed  (see  Fig.  3(b)].  There  is  no 
collapse  of  the  shear  modulus,  which  is  reduced  only  by 
an  amount  comparable  to  the  reduction  in  density. 

In  summary.  Just  above  the  self-diffusion  temperature 
and  well  below  the  melting  temperature,  only  a  very  low 
concentration  of  isolated  vacancies  can  be  sustained  in  a 
monolayer.  What  happens  to  the  vacancies  depends 
strongly  on  the  range  of  the  interparticle  interaction.  A 
short-range  interaction  favors  cavitation.  Monolayers 
with  long-range  interparticle  interactions  permit  the  heal¬ 
ing  of  the  monolayer,  leaving  behind  a  dislocation  dipole 
and  a  disordered  state.  The  speed  of  this  healing  is  re¬ 
markable,  of  the  order  of  nanoseconds  when  rare  gas  pa¬ 
rameters  are  used. 
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The  latter  mechanism  may  be  applicable  to  the  melting 
of  some  physisorbed  monolayers,  in  particular  those  with 
a  weak  substrate  corrugation,  such  as  Xe  on  Agil  I  iK  In 
an  experimental  environment  thermal  vacancies  are  un¬ 
avoidable.  They  may  permit  a  melting  transition  involv¬ 
ing  the  creation  of  dislocation  dipoles  and  their  dissocia¬ 
tion,  a  process  proposed  by  KTHNY  [2]  and  which  some 
experiments  support.  Existing  simulations  ignore  in- 
grown  vacancies  and  do  not  see  such  a  mechanism  |S,6). 
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Abstract:  The  methods  of  Boltzmann's  classical  statistical 
mechanics  have  been  used  to  calculate  the  hardening  at  OK  due  to  a 
random  dispersion  of  spherical  particles  with  finite  size. 


Introduction 

Previous  analytic  methods  for  treatment  of  dispersion 
hardening  have  been  confined  to  point  (i.e.  dimensionless) 
obstacles.  In  the  work  described  here,  these  methods  are  extended 
to  deal  with  the  case  where  the  particles  are  of  finite  size, 
spherical  and  of  radius  r. 

We  seek  to  determine  the  yield  strength  at  OK  of  a  material 
hardened  by  the  presence  of  spherical  particles  which  act  as  simple 


1 


blocks  to  the  motion  of  dislocations.  The  analysis  is  based  on  a 
consideration  of  the  stability  of  a  sequence  of  lengths  of  a  single 
dislocation  which  has  impinged  on  spherical  particles  while 
traversing  a  slip  plane. 


Analysis 

We  assume: 

that  when  acted  on  by  a  uniform  stress  o,  dislocations  with  Burgers 
vector  b  are  bent  into  arcs  of  circles  of  radius  R, 

(1) 

a 

where  p  is  the  shear  modulus  of  the  deforming  material; 
that  impenetrable  particles  are  randomly  distributed  and  are  of 
strength  F  which  can  be  expressed  in  terms  of  the  dislocation  line 
tension  ~  pb^/2  and  the  breaking  angle  a  (a)  (a  is  the  radius  of  the 
circle  of  intersection  between  particle  and  slip  plane)  by 

F  =  pb^sin  a  (a) 

where  each  breaking  angle  (a)  relates  to  an  individual  particle, 
being,  in  particular,  a  function  of  the  radius,  a,  of  the  circle 
of  intersection  of  the  operative  slip  plane  with  the  particle,  and; 

that  the  yield  stress  is  the  largest  stress  for  which  it  is 
possible  to  have  a  stable  dislocation  join  an  infinite  number  of 
particles  along  a  line  which  is  essentially  straight  and  hence  does 
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not  cross  itself. 


Towards  the  determination  of  this  stress  we  shall  use  the 
theoretical  analysis  first  introduced  by  Hanson  and  Morris  [1]  and 
subsequently  modified  and  employed  by  Labusch  [2]  and  by  Louat  (3) 
to  determine  the  yield  strength  due  to  a  dispersion  of  point 
obstacles  characterized  by  the  same  maximum  restraining  force: 

F  -  nij^sina^,-  (2) 

In  this  case  the  breaking  angle  a  is  the  same  for  all  obstacles. 

This  analysis  rested  on  the  rolling  circle  technique  employed 
by  Foreman  and  Makin  [4]  in  their  numerical  investigation  of  this 
hardening.  Thus,  it  was  argued  (see  fig.  la)  that  a  dislocation 
stable  at  a  sequence  of  positions  -  -  -,k-2,  k-1,  would  also  be 
stable  at  the  next  position  k  provided  the  particle  associated  with 
the  k+l’^*'  position  was  centered  anywhere  in  the  region,  as  shown 
hatched  (fig.  lb),  traversed  by  a  rolling  circle;  this  region  is 
bounded  by  an  angle  0  ^  2a.  Further,  for  a  concentration  of 

obstacles,  c,  the  probability  that  the  point,  k,  lies  in  a 
particular  area  AA  of  A  is 

provided  c  6A  is  sufficiently  small. 


Fig.  la.  Rolling  circle 


Fig.  l.b.  Relationship  between 
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sweeping  area  A. 


t„  and 


In  order  that  such  a  sequence  of  lengths  represent  a 
dislocation  line  it  is  necessary  that  they  not  cross  themselves. 

Here,  we  shall  adapt  the  treatment  of  Hanson  and  Morris  [1] 
to  the  case  where  the  strength  of  a  particular  obstacle  depends  on 
the  radius  of  the  circle  of  intersection  between  the  spherical 
particle  concerned  and  the  slip  plane.  To  proceed,  we  note  first 
from  fig.  (1  a.),  in  the  case  of  point  obstacles  that,  following 
the  path  of  a  rolling  circle,  the  angle  turned  between  successive 
tangents,  tn,  is  equal  to  the  angle,  <|>,  which  gives  the  direction 
of  the  trailing  circle  edge  at  the  point  of  rotation  (see  fig  (1 
b) ) .  It  is  then  readily  seen  that  in  the  case  of  an  obstacle  with 
radius  ,  a,  the  angle  turned  becomes: 

^  =  4>  (l  +  a/  R)sT<^, 
where  R  is  as  defined  in  (1) . 

Then,  if  the  number  of  points  be  N  and  if  the  angle  4>  occurs 
Ui  (i  =  l,n)  times  we  have; 

N  =  (3) 

i-l 

and  the  total  angle  turned  is 

S  ®  (4) 

i«l 

If  N  is  indefinitely  large,  but  9  finite,  the  average  value  of  <|>i 
is  zero.  In  this  case  the  whole  dislocation  line  is  either  straight 
or  contains  equal  numbers  of  complete  clockwise  and  anticlockwise 
rotations.  The  latter  configurations  are  situations  of  low 
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probability  which  we  can  neglect. 

The  procedure  for  the  determination  of  the  thermodynamic 
probability  is  akin  to  that  employed  to  determine  the  Maxwell- 
Boltzmann  distribution  of  kinetic  energies  in  an  ideal  monatomic 
gas.  Thus,  if  there  are  n  partitions  at  each  of  N  points  the 
thermodynamic  probability  that  the  i*"'  region  is  occupied  n^  times 
is : 

(N  !)  np/^ 

P=  - (5) 

n  n  ! 

i-1 

where  p^  is  the  occupational  probability  of  the  i’^*'  region  .  To 
determine  Pi  we  divide  the  swept  area.  A,  of  fig.  (lb)  into  n 
regions  of  area,  AJ'i  (i  =  l,n)  each  characterized  by  a  value  of  (J), 
say  and  see  that  when  the  concentration  of  second  phase 
particles  is  C  and  provided  the  area  AAi  is  sufficiently  small,  the 
probability  that  a  point  k  lies  in  a  particular  area  AA^  is: 

-3  n  f  SAj 

p_j  3  =  l-e  ~  CAA_f  (6) 

Here,  we  have  expressed  C  in  terms  of  the  volume  fraction,  f,  of 
second  phase.  Thus,  we  have  allowed  that  the  number,  M,  of 
particles  per  unit  volume  is  such  that 

~3 

where  f  is  the  volume  fraction  of  the  dispersed  phase.  It  then 
follows  that  the  concentration,  C,  of  circles  of  intersection  of 
particles  by  a  slip  plane  is: 
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C=2r  M  = 


3f 


2iir* 

On  the  basis  that  the  most  probable  configuration  of  a  system 
having  a  large  number  of  parts,  is  that  which  actually  occurs,  we 
see  that  to  determine  this  configuration,  we  have  only  to  maximize 
it  probability,  P,  subject  to  the  constraints  offered  by  (3) 
and (4) .  To  proceed,  we  take  advantage  of  the  facts  that  N  and  n^ 
are  both  large  and  employ  Stirling's  approximation  for  the 
factorials  in  (5)  and  then  set  the  variation  in  In  P  to  zero.  We 
find: 


6  Inp 

while  (3)  and (4)  become: 


a  o 

S  6n^  +  In  pj-  S  In 

i*l 


2  8nj=  0, 
i-1 

(7) 

n 

s  fixjj  4>i  =  0 . 

(8) 

Whence  introducing  the  Lagrangian  multipliers  k  and  y  we  have 

In  Pj-  In  nj+  Y  0 

so  that,  using  (3) 

.Hi 


n,.  = 


_  Np^e' 


i-1 


(9) 


Then  from  (4)  we  have 

Then,  on  substituting  for  n^  in  the  logarithm  of  P  and  after  some 
manipulation  we  find  that: 
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Nil  <|>i  Pi  =  ®  S, 

i-1 

where 

E4>ie^*  =  S.  (10) 

i-i 


P  =  s"e-^  (11) 

Invoking  (6)  and  (9)  and  proceeding  to  the  limit  of  small 
quantities  fiA^  ,  (11)  becomes: 

P  =  e'^lcj'e^dA]*'. 

We  have  now  to  determine  S  and  A,  .  Here,  where  integration  is  over 
the  area,  A,  indicated  in  fig.  (la)  we  have,  for  a/2R«l: 

2i2  r(2c-arcalni  /2S) 

S  =  cj  {1  +r  )  dl  J 

0  T(«rc8ln  i  /ZR) 


C  R^T  [e^**‘'^-l) 

X  (AD  2+1 


(12) 


The  determination  of  A.  is  assisted  by  the  realization  that 


dins  .  1  dS  .  ,13, 

dX  S  dX  S  N  ^ 


and  that  X  must  be  such  that 


dX 


=  0. 


Nevertheless,  this  question  can  only  be  resolved  numerically.  In 
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the  special  case,  unrealistic  here,  but  appropriate  for  point 
obstacles,  in  which  a  «  1  we  find,  independent  of  the  value  of  T, 
that 

X  =  1.41  /  a. 


In  comparison  with  the  cases  treated  previously  [1-3]  there 
are  additional  complications:  the  obstacles  have  finite  radius  and; 
as  we  shall  now  see,  have  variable  strengths  as  a  consequence  of 
having  variable  radii.  Bacon  ,Kocks  and  Scattergood  [4]  have  shown 
through  numerical  studies  that  the  Orowan  stress,  t,  for  a  linear 
sequence  of  impenetrable  circular  blocks  having  radii  a  and  inter¬ 
center  spacings  L  is  given  by: 


where. 


X 


(L-2a] 


\ib  since. 


_X  _  X 

X  ~  2a 


(14) 


In  this  expression,  we  have  determined  the  breaking  angle  in  the 
approximation  (2)  by  equating  the  Orowan  stress  to  that  found 
numerically  [5] .  Thus,  the  specific  strength  of  these  obstacles  is 
variable  and,  as  we  have  supposed,  can  be  denoted  by  the  term: 


sina 


ln(^ 

4« 


(15) 


This  supposition  may  be  justified  from  the  observation  [5]  that,  at 
breakaway,  an  obstructed  dislocation,  with  symmetrical  pinning,  has 
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swept  an  area  approximating  that  of  a  segment  of  a  circle  which 
subtends  an  angle  2a  at  its  center.  That  is  to  say,  taken  as  a 
whole,  the  dislocation  bends  in  much  the  same  way  as  in  the 
idealized  model  and  thus  into  a  circle  of  radius  R,  given  by  the 
relation: 

R  = 

t 

The  calculation  of  the  yield  strength  for  the  circumstances 
specified  here  is  complicated  by  three  features .  These  are:  the 
obstacles  are  not  distributed  along  a  straight  line;  they  are  not 
regularly  spaced  and;  their  strengths  are  variable.  The 
consequences  of  the  first  two  particulars  have  been  considered  in 
the  foregoing,  we  now  examine  the  effects  of  the  variation  in 
specific  strength. 

To  progress  we  suppose,  as  indicated  earlier,  that  the 
particles  are  spherical  of  radius  r  and  that  the  concentration  of 
the  particulate  phase  is  f.  We  then  find  that  the  mean  spacing 
between  particle  centers  in  a  slip  plane  is 


Again,  the'  mean  value  of  the  radii  of  intersection  of  such  a 
plane  with  the  particles  is  readily  found  to  be: 


Then  towards  the  evaluation  of  the  quantity  sin  a  as  defined 
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in  (15)  we  follow  Friedel  and  suppose  that  at  yield,  dislocations 
are  bent  into  arcs  of  circles  of  radius  R  such  that  the  areas  swept 
by  the  dislocation  in  circular  segments  between  obstacles  contain 
about  one  obstacle.  Taking  the  critical  angle  turned  by  the 
dislocation  at  an  obstacle  as  26,  the  area  of  a  swept  segment  is: 

4R20V3. 


Further, 


^  {L-2a)^Qf  ^  2ita^ 
3  3  3 


Thence  we  see  that 

L-2a  _  rZlT 
a  "  V  f  0 


so  that 

Thus,  the  mean  value  of  the  quantity  X/b  defined  in  (14)  is 


2a 

jb 

provided,  as  is  usually  the  case,  f0  <  »  0.1,  and  (15)  becomes: 


4lC 


=  sina. 


That  is  to  say,  in  the  circumstances  considered  here,  the 
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strength  of  a  particular  obstacle  of  radius  a  is  denoted  by 


a-  sin  a  =  — — (16) 

41C 

An  immediate  consequence  of  the  foregoing  is  that  the  radius  R  »a 
since  it  is  significantly  larger  than  L  unless  0  >  60“. 

We  are  now  in  a  position  to  consider  the  consequences  of  the 
simultaneous  existence  of  obstacles  of  differing  strengths. 
Reference  to  (12)  and  the  realization  that  the  analysis  which 
proceeds  it  is  essentially  linear  in  the  concentration  of 
obstacles,  shows  that  we  may  define  a  parameter  for  each  group 
of  obstacles  characterized  by  a  strength  and  a  concentration  Cj 
and  then  consider  the  quantity: 

Cy(e^*^^-1) 
o  =  _ 

®  X(X^T+1) 

which  in  the  limit  of  a  continuous  distribution  of  obstacle  sizes 
becomes : 


=  -TTrr— r  (e^-^^-Dda.  (17) 

X(X2r+i){ 

where  a„  is  the  maximum  strength  of  an  obstacle,  and  C,  is  the 
concentration  of  obstacles  of  strength  a. 

In  the  case  under  consideration,  where  the  particles  are 
spherical  and  all  of  the  same  radius,  r,  it  is  easily  shown  that 
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the  fraction,  f(a)  of  circles  of  intersection  of  having  radii  in 
the  range  a  to  a  +fia  is 


B[x,y) 


.ru)r(y) 

r(x+y) 


1  Z-1 

r(2)  ~  2«  ^  e'*  z  2  . 

Whence  we  find  that  (19)  is: 


j .  uliIul. 

(A+2ll)^ 

and  that  the  inequality  (18  )  becomes: 

X 

r2  c  r'^ 
2va^  x{x^+i) 

<So 

-K*— 

^^2  ^  ^  2i(^  2lC-ft 

2V^+2'it  X  (X*+l) 


Inspection  shows  that  these  bounds  are  nearly  coincident  so  that  we 
may  concentrate  our  attention  on  the  simpler  of  the  two. 

Towards  the  determination  of  X  we  have  as  indicated  earlier: 

dln5o  _  1  d5o  ,  ^  ^ 

dX  "  5o  dX  '  5o  N' 

a  quantity  whlcn  tends  to  0  as  N  tends  to  infinity.  Accordingly 
differentiating  with  regard  to  X  we  find  that  we  require 
and  so  from  (20) 
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=  0. 


d5o 

dX 


=  0 


1  _  1  _  2X  _  _1 

27C  b  X+2ii  X^+1  ^ 

In  the  range  of  interest  in  which 

102  s  i  10® 
b 

we  find  that 

2^X21. 

Finally  we  require  as  a  condition  for  yield  that  So  =  1 . 
Accordingly  setting 

R 

O 

in  (20  )  and  substituting  for  C  in  terms  of  f  we  have 

3\i^b^f{^)  2* 

o2  =  -  - -  (21) 

4ji  r2^X+27t  X(X2+i) 

Since  X  is  largest  when  r  is  least  we  find  that  there  is 
little  variation  in  the  term  (Zr/b)*''^*  with  variation  in  X.  Overall, 
for  a  constant  value  of  f,  the  yield  stress  varies  nearly  with  the 
inverse  of  the  particle  size.  It  may  be  noted  that  Freidel's  [6] 
conclusion  that  the  yield  stress  varies  with  the  square  root  of  the 
concentration  of  point  particles  is  retained  to  the  situation  where 
the  particles  have  finite  radius. 

It  may  also  be  noted  that  using  the  relation: 

(21)  may  be  written  in  the  more  direct  form: 
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r2  =  2nr^ 

^  3f  ' 

02  =  - £ -  . 

X(A2+i) 

Expressed  in  this  way,  we  see  that  the  yield  stress  varies 
essentially  with  the  inverse  of  the  mean  particle  spacing. 


Conclusions : 

An  analytical  treatment  which  proceeds  on  the  same  lines  as 
those  previously  employed  to  evaluate  the  hardening  effect  of  a 
dispersion  of  point  obstacles  has  been  used  to  deal  with  the  case 
where  the  particles  have  finite  radius.  Qualitatively,  the 
principal  conclusions  are  that  the  predicted  hardening  should 
increase  with  the  square  root  of  the  concentration  of  particulate 
phase  and,  for  a  fixed  concentration,  approximately  with  the 
inverse  of  the  particle  diameter.  Further,,  the  individual 
dependencies  of  strength  on  concentration  and  particle  size  can  be 
subsumed  into  a  single  parameter:  the  inverse  of  the  spacing 
between  particle  centers. 
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